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The null infinity limit of the gravitational energy-momentum and energy flux determined by 
the covariant Hamiltonian quasi-local expressions is evaluated using the NP spin coefficients. The 
reference contribution is considered by three different embedding approaches. All of them give the 
expected Bondi energy and energy flux. 
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I. INTRODUCTION 
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It is well-known that, as a consequence of the equivalence principle, gravitational energy cannot be localized (see 
e.g. Q §20.4). An alternative idea is quasi-local, namely quantities associated with a closed two-surface Q. During 
the recent decades there have been numerous intensive efforts made in the search for a better definition of quasi-local 
energy (as well as momentum and angular momentum) for gravitating systems, with the goal of obtaining quasi-local 
quantities which can provide a description of the gravitational field more elaborate than that given by the total 
quantities. A very nice review on the development and applications of quasi-local quantities in general relativity 
can be found in [3|. We are interested in testing certain quasi-local expressions for energy-momentum and energy 
■ flux obtained from the Hamiltonian boundary term using the covariant Hamiltonian formalism applied to gravity 

o: sua 

Some basic criteria are usually presumed for a physically reasonable definition of quasi-local gravitational energy 
0,13. One of the most important is to consider the asymptotic behavior of the quasi- local energy when the two surface 
approaches null and spatial infinity. The first aim of this work is to check, using the NP spin coefficient formalism 
0, the null infinity limit of the quasi- local energy-momentum associated with the boundary of a finite region for a 
5^ , certain covariant Hamiltonian quasi-local energy-momentum expression. 

We are interested in gravitational quantities like energy in the quasi-local sense, i.e. within a finite region of the 
space-time. There is an important issue concerning how much e nerg y flux flows into and out of the considered region. 
There continues to be considerable interest this topic 0, ll'2l ll.3|. In this work we also investigate a natural 
expression for the energy flux associated with the aforementioned quasi-local energy In order to test whether the 
definition is suitable, we again look to the the null infinity limit using the spin coefficient techniques. In this case, the 
gravitational energy flux is expected to be given by the well known Bondi energy loss formula. 

In order to have a reasonable definition of gravitational energy, the choice of reference plays a essential role. It is 
well known that the technique for choosing the reference is an important unsolved problem for the quasi-local energy 
issue. The ambiguity comes from how to embed the reference configuration into the physical space-time. Here we 
consider three different embeddings. 

This paper is organized as follows: In the next section, we start from a basic review on the asymptotic behavior 
of space-time near null infinity including a discussion of the expansion of the Newman-Penrose coefficients which is 
more complete than the well-known ones in |gL Il4l ITU Il6l Il7j . Our main result is contained in three parts: In Section 
III, we review the covariant Hamiltonian formalism and the associated boundary term approach to quasi-local energy 
momentum and energy flux. We then re-write the energy-momentum expression in terms of the Newman-Penrose 
formalism. In section IV, we find the asymptotic behavior at null infinity under three different embedding methods. 
From the results, we find that the Brown-Lau-York (BLY) embedding [l(| directly gives the standard Bondi mass 
aspect, while two other embeddings include an additional term which, however, vanishes upon integration. We identify 
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the source of this difference. In section V, we look at the direct definition of the energy flux, and consider its null 
infinity limit to test whether this definition is reasonable. The detailed calculation is divided into several subsections. 
In the first subsection, we calculate the purely physical part of the energy flux and, in later subsections, we will 
consider the three different embedding methods which have been mentioned. All three types of embedding give the 
standard Bondi energy flux. In section VI we test, using the spin coefficient formalism in the null infinity limit, the 
new Hamiltonian identity based expression for energy flux. It directly yields the expected Bondi energy-flux value. 
Section VII includes our concluding discussion. 



II. ASYMPTOTIC BEHAVIOR OF SPACE-TIME NEAR NULL INFINITY 



We first review the asymptotic behavior of space-time near null infinity. Here we are interested in the cases in 
which the space-time is asymptotically flat. There have been many intensive investigations of this subject in the past 
decades. We will follow the method initiated by Newman, Penrose, and Tod [l5lll6| . Many additional useful results 
have been discussed in [H 03, El E E3 

In accord with Penrose's conformal compactification method, i.e. the Penrose diagram, we assume that (M, g^, f2) 
is the conformal compactified manifold of a physical space-time (A/, g p „) via a conformal transformation g^v — g^v 
Suppose S is a section of future null infinity X + , we choose Bondi coordinates near X + in the following way: (#, (j)) are 
spherical coordinates on 5* and u is the affine parameter generating I + such that u — on S. If is a null vector onX + 
which is not tangential, namely fc M ^ T(I + ), then the null geodesies generated by fc^ will take the coordinates (u, 9, </>) 
into the physical space-time (M, g^ iv ). Moreover, on each null geodesic ip( r ) — { u = const, 9 = const, <\> = const}, its 
affine parameter r, in the sense of the physical metric g^u, can serve as the forth coordinate (in the neighborhood of 
X + , we assume the null geodesies are complete). Thus we have constructed coordinates (it, r, 9, <j)) for the physical 
space-time (they are like Bondi's coordinates except that r is an affine parameter rather than a luminosity distance, 
a small difference asymptotically). 

In addition to this coordinate construction, we impose, as has been introduced in |l9U20j. certain null frame gauges 
choices; stated in terms of the spin coefficients they are 

p — p — fi — p,~K^e = r — a — /3 = 7f — a — f3 — 0. (1) 

The fact that p and \i are real functions is insured by the null frame m and m being always tangent to the two sphere 
r = const and u — const. The meaning of the other gauge choices is that I is a null geodesic and the Lie transport of 
both n and m along / are tangent to the sphere. 

We use the standard notation for derivatives, D := l a S7 a , D' :— n a W a , S := m a W a , 5 '■= m°V a , and the usual 
complex coordinate on the two sphere, C : = cot | e 1 ^. The null tetrad can be chosen as 

1 - I W 
d d d-d 

" = 1Tu +u T» + x 1k +x TC (3) 

- = £^> ( 4 ) 

where U (real) and X, £,77 (complex) are undetermined functions. 

For simplicity, in this paper we will only focus on the vacuum case. Due to the asymptotic flatness, the behavior 
of the Weyl curvature satisfies the "peeling off theorem" 0,^3; i- e - 

*,~0(r 1 - 5 ), 1 = 0,1,2,3,4. (5) 

Near null infinity, we expand all quantities in Taylor series with respect to 1/r. Using the Newman-Penrose (NP) 
equations, the asymptotic behavior of the spin-coefficients are [T^.l20| 

1 Icr°| 2 

+ 0(r- 5 ), (6) 
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■ 1(71(7 -^ + 0(r-% (7) 
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f 3 = -—- — + °( r )> (9) 
T = n = ^ + 0(r~ 3 ), (10) 

where a°(u, 6, <j>) and \&o(u, 0, </>), the leading order terms of a (of order r~ 2 ) and (of order r -5 ), are free functions. 
Moreover the variable a is an abbreviation for a := ^/jC the spin- weight operator ft is defined by ^ := + 

2sa° acting on a variable with spin- weight s 0, and P := 1 + CC = sin -2 |. For calculating the energy flux, we 
also need the asymptotic expansion for the other NP coefficients. From the vacuum NP equations, after imposing the 
gauge conditions (JTJ, we have 

D-f = (t + n)a + (f + 7r)/3 + tit + #2- (11) 
Using the results in Eqs. (|BT llO[) and the asymptotic condition (JSJ), it is straightforward to get 



7=^ + 0(r- 3 )= ~ ™" ' r V U " ^ + 0(r- 3 ). (12) 

In order to get the asymptotic expansion of the remaining NP coefficients we need to have some more control on 
the null tetrad. From the commutation relations, one can easily derive the null tetrad control equations 

DU = -(7 + 7). (13) 

DX = (f + tt)£ + (r + W)rj, (14) 

D£ = ti + cm, (15) 

D V = vZ + pq, (16) 

which lead to the following asymptotic behavior of the null tetrad 

U = -l + '-^ + Oir-*), (17) 
2 r 

X = -^ff^ + 0(r- 3 ), (18) 

£ = 7 + ^^ + 0(r- 5 ), (19) 

V = -^f + 0(r- 4 ), (20) 

where the numerical factor — | in (|17fl and the value of £° (£° := -^) are specified by the result from Minkowski 
space-time. The asymptotic behavior of p and A can now be retrieved from the following NP equations 

Dp-Sir = (pfJ, + aX) + |tt| 2 - (a-(3)n + ^ 2 , (21) 

DX-5n = (pX + ap) + 7T 2 + {a - (3)n, (22) 

5t - D'a = {pa + Xp) + t 2 - (a - /3)t - (3 7 - 7 )cr. (23) 

The last equation is needed to determine the value of the leading order term of A. Using the obtained results (|17H20|) 
and the property that the spin weight of % is —1, it is straightforward to get 

* = 4- * S+ ^ g W °). m 

A ^ + i£!_p£! +0(r - 3) . (25) 

where the dot means derivative with respect to u and the numerical factor — ^ in the leading term of p is specified 
by the Minkowski space-time result. (It will turn out that p and A will play the key roles in our results.) 
Finally the expansion for the coefficient v can be obtained from the NP equation 

Dv - D'tt = (tt + f)p +(n + t)X + (7 - j)ir + * 3 , (26) 
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which gives 

^ + *W-). (27) 



However, the leading term of the NP equation 

8X-5fjt = {a + (3)n+{a-3f3)\-y 3 , (28) 

shows that "fg = —0 o a° (the spin weight of a is —2). Therefore the leading order of v indeed is r~ 2 which does not 
make a contribution in the later calculation. 

Moreover, for the later calculation, we still need the asymptotic properties of the induced volume element (2- 
dimensional) on the considered sphere S 2 . The induced metric ^ds 2 on the sphere S = {u = const, t = const} 
asymptotically should be 

Wds 2 = r 2 (d6 2 + sw 2 9d4> 2 )+0(r). (29) 

The induced volume element on S 2 can be obtained from the volume element e := ilAnAmA m, in our case 
e = ^J~^gdu A dr A dd A chf», by 

2 _ (d_\ a (d_ S " 

^cd — ^abed I ^ I I ^. 

\OU J \or 

= i m A m + 0{r) 

= r 2 (l- sin 0d9Adcj) + 0( r - 2 ). (30) 

In the derivation we have used the relation = — sin 2 | e~^(£g + -^g-§^)- 



ILL QUASI-LOCAL ENERGY-MOMENTUM AND ITS NULL INFINITY LIMIT 

There have been many proposals regarding quasi-local quantities. It should be noted that there is as yet no consensus 
regarding what approach should be used or even what are the proper criteria [3j. However it has been argued that 
the Hamiltonian approach, which has been used by many researchers and which we adopt here, has certain merits 
(see e.g. 0,la,Sllfllllll2H23jl23|)- For a general region S (finite or infinite) the Hamiltonian, 

H(N, E) = / N a H a + I B(N), (31) 

which displaces the region along a vector field N, includes not only an integral of a density over the 3-dimensional 
region but also an integral over its closed 2-surface boundary S = 9E. For Einstein's general relativity (as well as 
any other geometric gravity theory) the Hamiltonian densities H. a are proportional to certain field equations — the 
initial value constraints — and so the value of the Hamiltonian, E(N, S), is determined purely by the integral of the 
boundary term. For appropriate choices of the displacement N on the boundary, this Hamiltonian boundary term, for 
any gravitating system, determines the quasi-local values: in particular from a suitable time-like translation, the quasi- 
local energy and from a suitable spacclikc translation, the quasi-local linear momentum. The approach is quite general; 
it can incorporate in particular not only all the Noether charge expressions but also all the traditional pseudotensor 
expressions (and thereby it rehabilitates this often discredited approach) while taming the notorious ambiguities: the 
choice of boundary expression is linked, via the boundary term in the variation of the Hamiltonian, with the choice of 
boundary condition, while the reference frame ambiguity can be associated with a choice of boundary reference values, 
which determine the choice of vacuum or ground state for the system. In this way the traditional ambiguities can 
be given a clear physical and geometric significance ||| l2ll l22j . Within the covariant Hamiltonian formalism certain 
covariant symplectic expressions for the conserved gravitational quantities have been proposed [4J, |a, ||J . When we look 
at GR from this perspective one of these expressions stands out as being suitable for most applications, (among other 
virtues it has an associated positive energy proof Here we consider only this particular expression, specifically, 

in geometric units 

E(N, S) = ^i- £ (Ac^ b A i N r, ab + V^At^) , (32) 
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where Auj ab := uj ab — 10 ab is the difference between the orthonormal frame connection one-forms (i.e., the Ricci rotation 

o o 

one- forms) and their reference values, V and 7V a are the connection and the displacement vector in the reference space- 
time, rj ab := (l/2)e abcd '& c Atf d , A?7 ah := (l/2)e abcd {tf c A -d d - § c A § d ), with d a and § a being, respectively, the dynamic 
and reference orthonormal co-frames. 

The physical and geometric significance of this particular choice of Hamiltonian boundary term expression is revealed 
by the resultant boundary term in the variation of the Hamiltonian: 

6H(N, E) = / (field equation terms) - — <p i N (Auj ab A <577 a6 ). (33) 
Js lb7T Js 

This indicates that we should hold fixed on the boundary S the pullback of r] abl i.e., certain projected components of 
the coframe — thus, effectively, certain projected components of the metric (arguably the most natural choice). 

In the works already cited this quasilocal expression has been tested in various ways. Here we are concerned with 
the requirement that the value of the quasi-local energy-momentum and the energy-flux have the correct limit at 
null infinity. To this end we take S = <9£, the closed boundary of the three dimensional space- like region S, to be 
a two sphere which approaches in the limit null infinity. The Hamiltonian formalism with a boundary approaching 
null infinity has been considered from several perspectives, see e.g. ^JH^- A nice detailed discussion of the topic, 
addressing all of the important issues has been given recently 0] . 

By choosing the vector N, one can derive the 10 conserved quasi-local quantities for gravity based on the Poincare 
symmetry. In particular the covariant quasi-local energy and momentum associated with the time and space transla- 
tion asymptotic symmetries are 

PV = j s A lN " ^ + ^^Vab) ■ (34) 

Here the value of v labels what quantities are evaluated, for energy and 1,2,3 for the three components of momentum. 
In the asymptotically flat case, N v should be the translation part of the asymptotic Killing vectors. The translation 
part of the Bondi-Metzner-Sachs(BMS) group is well-defined; its expansion, in the leading order, is of the form 
N v = Nl a) + 0(r- x ), with v = 0, k = 1, 2, 3: 

^ = l = fo(*+h)+0{r-i), 



N i 0) = fk[^:-^z)=fk[n-^l)+0(r- 1 ), (35) 



du '"\ 2 

d d\ ( i 

where /„ = (1, — sin#cos</>, — sin#sin0, — cos#). 

The energy-momentum expression (|34H includes two parts which will be considered separately: the purely physical 
part and the part including the reference, i.e. 

P»=pZ hy +pf, (36) 

where 



Pl hy ■= J s (u ab A iN^ab) , pf := — J s (-> 6 A i NvVab + V^At^) . (37) 

We first evaluate the physical part and leave the reference part and the final results to the next section. It is 
important to keep in mind that the integral is evaluated on a two sphere with constant u and r, therefore the only 
contributing term is mAm. The gauge guarantees that the vector m is always tangent to the 2-sphere. Hereafter, 
we will only present the coefficient of the 2-form mAm, denoting this as "=" . 

From the expansion of the vector N, l|35|) . we realize that the the significant contribution for N = f v (n ± | /) is 

uj ab A i Nv Vab = fv {-«2(7 - 7) 1 A n ± i(a - (3 + f ± 2v) 1 A m =F i(a - + r ± 2D) 1 A m 

-i2(a - /3-?r)nAm + i2(a - /3-7f)nAm + i [2(fi + p) ± (p + p)] mAm}. (38) 

The imaginary unit i comes from the volume element, i.e. e = ilAnAmAm |2fij . However, if we only focus on the 
value on the two sphere boundary, the result is 

w ab A i^rjab S i.U [2( t i + p) ± (p + p)} m A m. (39) 
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IV. GRAVITATIONAL ENERGY-MOMENTUM IN DIFFERENT EMBEDDINGS 



Reference configurations play a crucial role in the expressions for gravitational energy and its energy flux (indeed 
for all the quasi-local quantities). There are two essential related issues: (i) a "suitable" reference configuration choice 
and (ii) a proper embedding into the physical space-time. 

O i 

There are two terms in the reference part. The NP formulation for the first one, u> a0 A iNVab can be easily read out 
from (|38(l by replacing all NP coefficients and the frames within the connection with their reference values. 

For an asymptotically flat space-time the choice of reference configuration is more or less unambiguous — the 
Minkowski space-time. In Eddington-Finkelstein coordinates (u, r, (9, <fi) , which are related to the standard coordinates 
by u — t — r, the first fundamental form of the Minkowski space-time is 



ds 2 



-du 2 - 2dudr + r 2 (d6 2 + : 



(40) 



It is easy to verify that the coordinates are just the Bondi coordinates. 

However, the issue of embedding the two surface S of physical space-time into the Minkowski space-time is not fully 
transparent yet. Various proposals have been made in earlier works. In this subsection, we consider several types of 
embedding used for calculating the reference part of the energy; later we will also use these embeddings for energy 
flux. We will see that the embeddings we consider actually give the same result for the gravitational energy flux at 
future null infinity, but not quite the same expression for energy itself. 



A. Holonomic embedding 



The simplest embedding technique is to identify the Bondi-coordinates with the Minkowski coordinates, i.e. the 
embedded surface So is just the standard coordinate round sphere. In this approach, the embedded null frame is 



8_ 

or 



d_ _ ld_ 

du 2 dr ' 

e _i * ( d id 



P d 



V2r \d0 sm6d4>) y/2r d(' 
Therefore the non-vanishing N-P reference coefficients are 

1 p 1 o 1 C 

P= — , f 1 = , a = -B= — 



(41) 
(42) 

(43) 
(44) 



The embedding of the time direction N is Nu = Ii/\ n + 1 )l) which is a Killing vector of the Minkowski space-time 



o , o 



(45) 



Properly we should calculate ui =T c , d c - However since the difference of the coframes -d c and $ c is o(r ) we can 
make an approximation and take for the reference part 

w ab M N riab = ifv 2(^ + fL)±(p + h mAm. 

o / D 0\ 

Moreover, it is straightforward to check, for Nv= fv \ n ±2 ^ ) + 0( r ~ )> tr±at 

= (d a fu)(n ±^ °l) b Ariab — ifv [(7 + ^)A(m A m)+ v A(l A m)— V A(l A m) iO^ 1 ) 



- Oir- 1 ). 

Finally, from the results l|39l 1451 146|) , the energy- momentum in the holonomic embedding is 

1 

i m A m 



(46) 



limp,, = lim— / f v 
1+ 1+ 16tt J s 



2(n + p- H - fi)±(p + p- P - p) 

= __L f Re(^ 2 + a ^ + f a°)f„dn 2 , 
** s 



(47) 
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where dfl 2 = sin6d6d(j). The integrand differs slightly from the usual formula for the Bondi energy-momentum; 
however the integral of the extra term vanishes at least for the energy — which is our real interest here — simply 
because $a° is of spin weight -1. Note that entirely analogous terms show up in equivalent calculations done directly 
in the Bondi-Sachs metric 0, |2g • 



B. 6 Murchadha-Szabados-Tod embedding 



In Ref. [23, 6 Murchadha, Szabados and Tod (OST) introduced another kind of embedding method. Let us 
consider the space-like region So m the physical space-time, 9So = S. We suppose that S is a topological two sphere 
and that isothermal coordinates globally exist on S. In these coordinates, the induced metric ^ds 2 on S is 



< 2 W = uj 2 {e'^'){de' 2 



,2 a* AA,'1 



), 



(48) 



where oj(9' , 4>') is a positive function on S. This two surface is embedded isometrically into the physical space-time 
(in Bondi coordinates) with 



u — const, 6 — 6', <p — 4>' , r = u>(6,(p). 
Based on this embedding, we define a coordinate transformation in Minkowski space-time, 

u^U, r^R + uj, 6^6, <t>^4>- 
Then the Minkowski metric becomes 



9 = 








I -1 -1 -d e cu 
-10 

-d e uj {R + uj) 2 
\ -d,puj (R + lu) 2 sm z 6 J 



The NP reference tetrad is chosen to be 

_ d_ 

1 ~ dR' 



d_ 

dU 



in 



-i<p 



|<W| 5 



d 



dgu d 



d<f,uj 



{R + u;) 2 J OR (R + uj) 2 d6 (R + u;) 2 sm 2 



V2(R- 



d_ 

06 



d 



sin 6 1 



where 6q = Sz + -rt-z Si ■ The associated dual reference tetrad is 



08 1 sin 8 dip 1 



1 = -dU, 
1 , 

n = — 1 
2 



\6ou\' 



(R + oj) 2 



dU - dR, 



-dU + — =e' 
V2(R + uo) V2 



-i(j> 



(R + Lu)(dd + ismdd(j)). 



From the above metric, a direct calculation gives the non-vanishing NP reference coefficients as 
1 



P 



R + u}' 



1 



\5 uj\ 



dnOJ + cot 9 deuj 



sin 2 8 



R + LU (R + lu) 3 



{R + uf 



2V2 
e^ cot § 
2V5R + lo'' 



cot I 25 uj 
R + uj + {R + u) 2 



(49) 
(50) 

(51) 



(52) 



(53) 

(54) 
(55) 

(56) 

(57) 
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7 = — 



\/2(i? + W ) 2 ' 
|<5 w| 2 icot^d<pui 



{R + lu) 3 (i? + cj) 2 sin( 



(58) 
(59) 



2{R + uo) 2 

e 2i<f> 

2{R + lo) 2 



2icot6d 4 ,u 2i 
. 2 - cotfldgw H — — r-^oeo^w 



+ w 



(5 2 cj — cot S uj 



R + uj 



V2(R + uf 



\SpUj\ 2 5 UJ 

R + uj 



— 5 uj Sodeiv + 



i cot Qd^ui \ id^u 



dgto + cot 9 doiu + — 



1 



sin' 6 



■ dluj 



(60) 
(61) 



The intrinsic geometry of the two sphere is preserved in the OST embedding, therefore A^j, = 0. Moreover, on 
the considered two sphere S, R = by the definition of the embedding and, from Eqs. and (pIU)) . we have 



T 0|2 



2r 



+ 0{r~ 2 ). 



(62) 



Therefore, this reference contribution differs from the result of the holonomic embedding only by higher orders of 1/r. 
Hence the energy-momentum at null infinity in the OST embedding is the same as in the holonomic embedding, i.e. 



lim n„ 



1 

47T 



Re (*g + a°a° + f a°) f v dfl 



(63) 



Again the result differs from the usual Bondi integrand by the same extra term — which makes a vanishing contribution 
to the energy when integrated over the 2-sphere. 



C. Brown-Lau-York embedding 

In the above two subsections, the two methods used both embedded the surface S into a standard light cone in 
Minkowski space-time, i.e. the light cone N is the light cone from one point. Brown, Lau and York gave another 
way to do the embedding near null infinity. This method considers a more general light cone. 

Suppose we choose a Bondi coordinate system (u,R,6,<f>) in Minkowski space-time. The asymptotic shear is <J°. 
We can also do the formal Taylor extension near null infinity as in Eqs. lj6*l ll(J|) . The only difference is that we have 

o 

i$fi= 0, n — 0, 1, 2, 3, 4. From the NP equations we have 

¥V=-0 o dJ°, #f = -c^°. (64) 

Because the spin- weight of d u a° is non-zero, the above results tell us that d u a° = 01- Insert this results into Eqs. 
(|6T llO|l . the N-P quantities in Minkowski space-time are then 



+ 0(i?- 5 ), 



1 k"l 

R R 3 

! ° 

0(R~ 5 ), 



R 2 R 4 

a la 

+ 0(R- 3 ), 



a a a + @ a 

a = 1 

R R 2 

a a a 

"-- 1>- — + "'^ 
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1 Ko 



2^0 



u i /nCD-3 



2R R 2 



0{R- 



v = 0{R~ 3 ). (65) 



The tetrad part is 



u = -i-2i±Ii + (fi-*), 

fl = --^+0(iT 4 ). (66) 
The intrinsic geometry of the 2-sphere is preserved in the embedding, i.e. 

TL = n, (67) 

where the two dimensional Ricci scalar TZ is given by 1Z = — 2p/i — 2pp + 2<rX + 2a X + 2^ 2 + 2^2- Consequently, this 
leads to 

- up + Xa + * 2 + *2 + Xa - ftp = ~ MP - Jip + °\<r + \e ■ (68) 
Using the Taylor expansion in Eqs. l(6"l ll0fl and Eqs. JBSJ, we find that the relation between the parameter r and R is 

R = r + k + 0{r- 1 ), with k = f a° + f £° - f Q cr° - f a° . (69) 

We choose cr°|s = cr°|s, where S is the section on I + and Sq is its image under the embedding. We have 

R = r + 0(r- 1 ) = r[l + 0{r- 2 )]. (70) 

As for the OST embedding, the two sphere geometry is preserved, therefore Ar) a i, = 0. Finally the gravitational 
energy-momentum in the BLY embedding is 

]jmp v = -^- [ Re(^° 2 + a°a ) f„dfl 2 . (71) 
j+ 4tt J s 

It is worth noting that the BLY embedding is a little neater, in that it directly gives the standard Bondi energy- 
momentum integrand, whereas there is an additional term in the other two embeddings (which vanishes upon integra- 
tion). That term is associated with the embedding methods in which the section S is not embedded into a standard 
light cone, generated by the null geodesies starting from a single point. This difference again shows us that keeping the 
inner geometry unchanged under the embedding is not enough to ensure physically reasonable quasi-local quantities; 
generally, as is especially clear from |27|. we need more restrictions. 

V. THE ENERGY FLUX AT NULL INFINITY VIA THE DIRECT METHOD 

In this section, we directly calculate the energy flux through a two sphere. For simplicity we choose the time-like 
translation to be N — n + ^l. Asymptotically, this agrees with the natural choice, the time translation of the BMS 
group at null infinity d u . The difference between the two vectors is asymptotically negligible. 

We consider the energy E = H(N,S). Suppose So is the space-like region that we want to consider, <9£ = S, 
and Ea( is the time evolution of So- The energy within the region during the time interval changes by the amount 
E(EAt) — E(Eq), hence a natural direct definition of the rate of energy change is 

E:= lim S(SA t )-£(So) 

At^O At v ' 
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Looking to the value of the Hamiltonian, taking into account the vanishing of the initial value constraints, from 1321) 
we straightforwardly get the quasilocal energy flux relation 



E = H{N, S) := ^i- £ £ N (Aw Qb A i N T) ab - V a N b Ari ab 



16tt 



£ N (uo ab A i N r] ab ) - £ N [uj ab A - £ N [V a N b ATj ab 



(73) 



The right hand side defines an energy flux expression F which includes two parts that will be considered separately: 
the purely physical part and the other part, which includes the reference, i.e. F = -F p hy + F le f. The reason for making 
such a separation is that the part including the reference, as we have already seen, depends on the embedding of the 
reference configuration into physical space-time. 
We first evaluate the physical part of the flux, 



£ N {io ab A i N riab) = \ £ i{u ab A i N r) ab ) + £ n {uJ ab A i N r) ab ). 



(74) 



We see that it is necessary to know the Lie derivative of all two form elements. However, we are only interested in the 
final results that can contribute to the integral: the terms proportional to m A m. After a straightforward verification 
the contributing terms are 



£i(m A m) = — (p + p) m A m, £ n (m A m) = (p + p) m A m, £ n (n A m) = — v m A m. 

Therefore, following the result of uj ab A ijv^ah m Eq. (|38|l . the first term is 

£i(to ab A i N r) ab ) = iD [(p + p) + 2(/x + p.)] m A m + i [{p + p) + 2(p + p)} £i(m A m) 
= * {D [{p + p) + 2{p + p)} -{p + p) [{p + p) + 2{p + /!)]} m A m 
= i [2Sir + 2Sn + TZ~ 2p 2 + 2\a\ 2 + 4|tt| 2 - 2(a - (3)tt - 2{a - (3)tt] m A m, 

where 1Z is the two dimensional Ricci scalar on the enclosed surface. Similarly the second term is 

£ n (uj ab A i N r/ ab ) ~ iD' [{p + p) + 2{p + fx)] m A m + i [(p + p) + 2(p + p)} £ n (m A m) 
-i2(a - ft - ir)£ n {n A m) + i2(a - /3 - ir)£ n (n A m) 
= * {D' [{p + p) + 2(p + fi)] + {p + p) {{p + p) + 2{p + p)} 
+2(a - (3 - n)& -2(a- - k)v) m A m 

5t + 5t + 28v + 269 - -TZ+ Ap 2 - 4|A| 2 + 2(7 + j)(p - 2p) + 2\ir\ 2 + (a- (3)w 



(75) 



(76) 



+ (a - /3)vf + 2v(a + 3/3) + 2v(a + 3/3) + 2(a - /3 - iz)v - 2{a - /3 - tt)v 



m A m, 



(77) 



Finally, the purely physical contribution to the energy flux is 
1 

"'phy 



— / £ N {uj ab Ai N r] ab ) 



I 

8^ 



Sn + 5n + 5v + 5v - ^p 2 + ^\a\ 2 



2p 2 + 2\ir\ 2 - 2\\\ 2 



+ (7 + j)(p - 2p) + v(a + 3/3) + D(a + 3/3) + (a - - ir)v - (a - /3 - tt)v 



m A m. 



(78) 



From the asymptotic expansion of all the NP coefficients l|6ll(JI 1121 1241 1251 1271) we find that all except one of the 
terms fall off as 0(l/r 3 ) or faster, only the — 2|A| 2 term contributes asymptotically. Hence the null infinity limit of 
the purely physical flux is 



lim F phy 



lim — , 

1+ 877 Js 

1 

4-7T 



r 

0|2 i 2 



0(r~ 



i m A m 



(79) 
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A. Holonomic embedding 



A straightforward calculation gives 



£ 



N 



(u ab M N Vab) = £N 



(p + p +2 p +2 p) i m A m + (a - 0) i n A m - (a - (3) i n A m + 0(r 3 ) 

i m A m, 



(p + p+2°H+2p){~p- ~p + p + p) + 0(r 3 ) 



A 



= 0. 

Therefore, the reference part in the holonomic embedding is 



loo _°°1° _°^1° °o 

-(7 + 7) i (m A m- m A m) + - V i (I A m— 1 A m) - - f i (1 A m- 1 A m) 



(80) 



Ft. 



1 



16tt Js 
Finally, we have 



{p+°p+2?i+2 p J ){--p- -p + n + p,) + 0(r 3 ) 



imAm= / Otr 3 )imAm. (81) 

167T J S 



lim F 

x+ 



47T 



|a°| 2 ^ 2 . 



(82) 



This is just the standard expression of the Bondi energy flux at X + |il IbH ITil l28| . 



B. 6 Murchadha-Szabados-Tod embedding 

The reference part in the OST embedding is 



£ 



N 



(> 6 



A iNVab 



o o 00- 0000 n -\- n 

(- 7f + a - + (- vr + a - + (p + p + p + ii)[-tL—H + M + ^) 



i m A m, 



( V^A^ ^ V a N h £N (Ar? a6 ) ^ - (7 + 7)£jvA(m A m)+ v £jvA(1 A m)— v £^A(1 A m) 



(83) 



The image of S is the two sphere in Minkowski space-time such that U — const, and R = 0. From the above 
calculation, we get the same result as before: 



i^ref = / 0{r~ 3 )i mAm. 
Js 

Consequently, for the energy flux we get the same result as Eq. 



(84) 



C. Brown-Lau-York embedding 

Following up on the result in the case of the energy-momentum calculation, the reference part of the energy flux is 

F Icl = / 0(r 3 )imAm. (85) 
Js 

Hence for the total energy flux we again we get the same result as Eq. I|82() . 



VI. THE ENERGY FLUX AT NULL INFINITY VIA AN IDENTITY 

In this section, we calculate the energy flux through a two sphere using an interesting formal Hamiltonian identity 
[29| derived in detail in 0- The identity is simply the analogue of the classical mechanics identity H = 0, which 
follows from SH = q k Spk — PkSqk by simply by replacing S — > d/dt. We can obtain it directly in essentially the same 
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way from (|33|l simply by substituting the time derivative operator £n for 5; then the field equation terms cancel 
identically (just as they did in the classical mechanics case) leaving 



E := H(N, Ti) = -<b t N (Atu ab A £ N i lab ) = d> (-i N Aw ab £ N r lab + Auj ab A £ N i N r lab ). (86) 
Js Js 

This is a general quasi-local formula for energy flux, applicable to the boundary of any region. In [3 it was tested in 
the null infinity limit using the Bondi-Sachs metric. 

Here we wish to transcribe this expression into the NP spin coefficient form and confirm that it gives the desired 

o , 

asymptotic results using that well developed technique. For our calculation here we take N — d u = n + ~l = 
n + hi + 0(1 /r), the reference geometry Killing field. 
Let us consider Eq. 1(86(1 term by term. The first term is 

i N Acj ab £ NVab = A^ b £ nriab + ^Auf£ lVab + ~Au? b £ n r] ab + ^A^ b £ lVab 

= 2A(7 + y)£ n (i m A m) - 2Av£ n (i m A 1) - 2Av£ n (i 1 A m) 
+2Af£ n (i n A m) + 2Ar£ n (i m A n) + 2A(j - j)£ n (i n A 1) 
+A(7 + i)£i(i m A m) - Av£ l (i m A 1) - Av£ l {i 1 A m) 
+A?£i{i n A m) + Ar£ ; (i m A n) + A (7 - j)£i{i n A 1) 
-A(e + e)£ n (i m A m) — An£ n (i m A 1) - A%£ n (i 1 A m) 
—An£ n {i n A m) — An£ n (i m A n) + A(e — e)£ n (i n A 1) 

+ ^[-A(e + e)£i{i m Am) - Air£i(i m A 1) - Air£i(i 1 A m) 

-An£i(i n A m) - An£i{i m A n) + A(e - e)£i(i n A 1)] . (87) 

Based on the asymptotic estimations given in the previous sections and the three considered embedding methods, all 
of these terms are of higher order than 0(X/r), hence they make no contribution to the flux asymptotically. Basically 
this comes about because (i) in our gauge e = K = 0, (ii) the Lie derivative introduces a factor of 1/r, (iii) the spin 
coefficients a, it, r, 7,1/ are 0(l/r 2 ), (iii) for p,a,/3,n the A operation removes their 0(l/r) part. 
Similarly, for the second term we find 

Auj ab A £ N iNVab = {-Auj^i rfl £ N i N 7] ab + Auj a i b i m £ N i N i] ab )m A m 

= 2[A(a + f3)irh£NiNVoi ~ (Ap)i rfl £ NiNVo2 - {A\)i r - n £ N i N r) 0S , 

+ (Ap)i fh £ N i N 7]i2 + (A(j)irn£NiNm3 + A (<5 ~ P)irn£ NiNll23 

-A(a + (3)i m £ N i N rioi + {AX)i m £ N i N r]o2 + (Ap,)i m £ N i N r) 03 
-(Aa)i m £ N i N rii2 - {Ap)i m £ N i N r] 13 + A(a - f3)i m £ jv^jv^s] m A m 

= -A\X\ 2 i mAm + 0(-). (88) 
r 

Again we used the asymptotic estimations and the three considered embedding methods. In this case we get a non- 
vanishing asymptotic contribution from A = 0(1 /r) 1)25(1 . Submitting these results into Eq. 1)86(1 . we find that the 
null infinity limit of the energy flux calculated from that Hamiltonian identity relation is 

limF = --!- / \&°\ 2 dn 2 , (89) 



X+ 47T 

which is, just as was found from the direct calculation 1(82(1 . the standard flux loss due to the Bondi news. 



VII. DISCUSSION 



We have tested certain expressions for the quasilocal energy-momentum and energy flux of gravitating systems. The 
expressions were obtained from the covariant Hamiltonian formalism. In this formalism the quasilocal quantities are 
determined by the value of the boundary term in the Hamiltonian. The variation of the Hamiltonian associates the 
choice of boundary term with specific boundary conditions. Thus the definition of the quasilocal energy-momentum of 
a gravitating system is linked to the choice of boundary conditions. The boundary term that corresponds to holding 
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certain projected components of the orthonormal frame fixed seems to be the best choice for most purposes. We have 
considered only that choice here (the values for certain other choices are given in Q). 

Energy flux can be computed in more than one way. On the one hand it can be obtained directly from the change 
in the energy expression. On the other hand one can use an interesting identity associated with the specific role of 
the Hamiltonian and its variation. Here we have evaluated the energy flux by both techniques. 

In strong field regions we do not have any sharp test as to what values we should find for energy-momentum and 
energy flux. Proposed expressions necessarily are first tested in the weak field linearized theory limits. Getting good 
values at spatial infinity is not the strongest test. The Bondi limit at future null infinity is more delicate. Here we 
tested our selected expression for energy- momentum and its associate energy flux in this limit. 

Technically we used a well known and well developed technique: the Newman-Penrose spin coefficients. We selected a 
suitable gauge and found that we needed certain quantities expanded in more detail than is usual |30| . In the quasilocal 
expressions it is necessary to select reference values which determine the "vacuum" or "ground state" . The natural 
choice is, of course, Minkowski space, but it is not so obvious how to embed the Minkowski space into the asymptotic 
part of the dynamic space. We considered three types of embeddings which have been used: holonomic, one due to 
O Murchadha, Szabados and Tod , and one due to Brown, York and Lau We found some interesting technical 
differences between the embeddings but in the end they all gave the same answer: namely the expected Bondi energy 
and the Bondi energy flux determined by the Bondi news. 

In the detailed calculation we noted that, at least in the selected gauge, the quasilocal energy was asymptotically 
determined by the deviation of the spin coefficient fi from its asymptotic Minkowski value and that the energy flux 
was determined by the spin coefficient A (in the notation of 01 t nese coefficients are —p' and —a', respectively). 

We have shown that the values of these expressions can be practically calculated in terms of the NP spin coefficient 
technique; the expressions were found to have the desired asymptotic values. Thus they satisfy an important criterion 
for quasi-local energy and energy flux expressions. 

Acknowledgements 

This work was supported by grants from the National Science Council of the Republic of China, X. Wu was 
supported by grant number NSC 92-2816-M-008-0004-6, C.-M. Chen by grant number NSC 93-2112-M-008-021, and 
J. M. Nester by grant number NSC 93-2112-M-008-001. 



[1] C. W. Misner, K. S. Thorne and J. A. Wheeler, Gravitation, (San Francisco: Freeman, 1973). 

[2] R. Penrose, "Quasilocal mass and angular momentum in general relativity", Proc. Roy. Soc. Lond. A 381, 53-63 (1982). 
[3] L. B. Szabados, "Quasilocal energy-momentum and angular momentum in GR: A review article", Living Rev. Relativity 

7, (2004), 4. [Online Article]: cited [1 April 2005] , |http:/ /www.livingreviews.org/lrr-2004-4| 
[4] C. -M. Chen, J. M. Nester and R. S. Tung, "Quasilocal energy-momentum for geometric gravity theories", Phys. Lett. A 

203, 5-11 (1995) [arXiv:gr-qc/9411048| . 
[5] C. -M. Chen and J. M. Nester, "Quasilocal quantities for GR and other gravity theories", Class. Quantum Grav. 16, 

1279-1304 (1999) larXiv^gr- qc/9809020]. 
[6] C. -M. Chen and J. M. Nester, "A symplectic Hamiltonian derivation of quasilocal energy-momentum for GR", Gravitation 

& Cosmology 6, 257-270 (2000) arXiv:gr-qc/0001088 . 
[7] C.-M. Chen, J. M. Nester and R. S. Tung, "Quasilocal energy flux", in preparation. 
[8] C. C. M. Liu and S. T. Yau, "Positivity of quasi-local mass II", arXiv:math.DG/0412292 . 
[9] R. Penrose and W. Rindler, Spinors and Space-time Vol. 1, 2, (Cambridge University Press, 1984, 1986). 
[10] J. D. Brown, S. R. Lau and J. W. York, "Energy of isolated systems at retarded times as the null limit of quasilocal 

energy", Phys. Rev. D 55, 1977-1984 (1997) arXiv:gr-qc/9609057 . 
[11] P. T. Chrusciel, J. Jezierski and J. Kijowski, Hamiltonian Field Theories in the Radiating Regime, (Lecture Notes in 

Physics, New Series M, Monographs vol M70) (Berlin: Springer, 2001). 
[12] J. W. Maluf, F. F. Faria and K. H. Costello-Branco, "The gravitational energy-momentum flux", Class. Quantum Grav. 

20 (2003) 4683 |arXiv:g r-qc/0307019 . 
[13] J. H. Yoon, "New Hamiltonian formalism and quasilocal conservation equations of general relativity", Phys. Rev. D 70, 

084037 (2004) |arXiv:gr -qc/0406047l. 
[14] J. Stewart, Advanced General Relativity, (Cambridge University Press, 1990). 

[15] E. T. Newman and R. Penrose, "An approach to gravitational radiation by a method of spin coefficients", J. Math. Phys. 
3 566-578 (1962). 

[16] E. T. Newman and K. P. Tod, "Asymptotically flat space-time", in General Relativity and Gravitation, Vol.2, ed. by 
A. Held (New York, Plenum, 1980). 



13 



[17] R. Geroch, A. Held and R. Penrose, "A space-time calculus based on pairs of null directions", J. Math. Phys. 14 874-881 
(1973). 

[18] M. Ludvigsen and J. A. G. Vickers, "Momentum, angular momentum and their quasi-local null surface extensions", J. 

Phys. A: Math. Gen 16 1155-1168 (1983). 
[19] R. M. Kelly, K. P. Tod and N. M. J. Woodhouse, "Quasi-local mass for small surfaces", Class. Quantum Grav. 3 1151-1167 

(1986). 

[20] A. Dougan, "Quasi-local mass for spheres", Class. Quant. Grav. 9 2461-2471 (1992). 

[21] C- C. Chang, J. M. Nester and C. -M. Chen, "Pseudotensors and quasilocal energy- momentum" , Phys. Rev. Lett. 83, 

1897-1901 (1999) arXiv:gr-qc/9809040 . 
[22] J. M. Nester, "General pseudotensors and quasilocal quantities", Class. Quantum Grav. 21, S261-S280 (2004). 
[23] J. D. Brown, S. R. Lau and J. W. York, "Canonical quasilocal energy and small spheres", Phys. Rev. D 59, 064028 (1999) 

arXiv:gr-qc/9810003\ 

[24] J. M. Nester, "Positive Energy Via the Teleparallel Hamiltonian" Int. J. Mod. Phys. A. 4, 1755 (1989); "A Positive 

Gravitational Energy Proof" Phys. Lett. 139 A, 112 (1989). 
[25] R. D. Hecht and J. M. Nester, "An evaluation of the mass and spin at null infinity for the PGT and GR gravity theories," 

Phys. Lett. A 217 81-89 (1996). 
[26] D. Kramer, H. Stephani, E. Herlt and M. MacCallum, Exact Solutions of Einstein's Field Equations, (Cambridge University 

Press, 1980). 

[27] N. Q Murchadha, L. B. Szabados and K. P. Tod, "Comment on 'Positivity of Quasilocal Mass' ", Phys. Rev. Lett. 92 
259001 (2004). 

[28] S. W. Hawking, "Gravitational radiation in an expanding universe", J. Math. Phys. 9 598 (1968). 

[29] The possibility of such an identity was realized long ago, see e.g. the remark on page 160 in J. M. Nester, "The Gravitational 
Hamiltonian", in Asymptotic Behavior of Mass and Space-Time Geometry, ed. by F. Flaherty, Lecture Notes in Physics 
202, (Springer, 1984) pp 155-63, but not followed up until recently. 

[30] Since completing our work, we have learned of another work with quite detailed spin coefficients: W. T. Shaw, Classical 
Quantum Gravity 3, 1069 (1986). 



14 



